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Through the exact solution of a two-electron singlet system interacting with a monochromatic 
laser we prove that all adiabatic density functionals within time-dependent density-functional theory 
are not able to discern between resonant and non-resonant (detuned) Rabi oscillations. This is 
rationalized in terms of a fictitious dynamical exchange-correlation (xc) detuning of the resonance 
while the laser is acting. The non-linear dynamics of the Kohn-Sham system shows the characteristic 
features of detuned Rabi oscillations even if the exact resonant frequency is used. We identify the 
source of this error in a contribution from the xc-functional to the set of non-linear equations that 
describes the electron dynamics in an effective two-level system. The constraint of preventing the 
detuning introduces a new strong condition to be satisfied by approximate xc-functionals. 



I. INTRODUCTION 

Despite of the success of linear-response schemes to de- 
scribe excitations of many electron systems, many physi- 
cal processes stemming from the interaction of light with 
matter are non-linear in nature (i.e. photo-reactivity, 
population quantum control, etc). A prototype effect, 
ubiquitous in many fields of physics, are Rabi oscilla- 
tions between the ground and an excited state when a 
monochromatic laser with a frequency close to the res- 
onance is applied. The dipolc moment of the system 
then shows fast oscillations with the frequency of the 
applied laser and an envelope that oscillates with the 
frequency attributed to the transition from the ground 
to the excited state (Rabi frequency). It is one of the 
few analytically solvable cases of non-linear electron dy- 
namics where the population of states changes dramat- 
ically in time. Thus, Rabi dynamics provides a firm 
ground to derive fundamental properties of the exact 
exchange-correlation (xc) functional of time-dependent 
density-functional theory (TDDFT) [1, 2]. This is one of 
the scopes of the present paper together with providing a 
sound description of laser induced population processes 
within TDDFT. We remark that it was recently con- 
cluded that Rabi oscillations cannot be described within 
a finite- level model in TDDFT [3] . The authors observed 
that the dipole moment in an adiabatic exact exchange 
(EXX) calculation shows the characteristic Rabi oscilla- 
tions while the density does not show a transition to the 
excited state. 

Here, we show that for all adiabatic functionals the 
observed oscillations in the dipole moment can indeed 
be described as Rabi oscillations, however, with a non- 
resonant laser frequency. The employed adiabatic ap- 
proximation results in a detuning, i.e. the system is 
driven out of resonance by the change in the Kohn-Sham 
(KS) potential which is due to the change in the density. 



The failure of adiabatic density functionals to correctly 
describe Rabi oscillations has some resemblance to well- 
known problems of describing double and charge transfer 
excitations [4, 5] and ionization processes [6] that were all 
traced back to the lack of memory in the functional. In 
the present work we use the example of Rabi oscillations 
to identify the physics behind the failure of adiabatic 
constructions to describe resonant non-linear dynamics. 
Our results show that the problems are not limited to 
adiabatic TDDFT but affect any mean-field theory. 

The paper is organized as follows. In Section II we de- 
rive the analytic solution for Rabi oscillations in a many- 
body interacting system and apply this solution to com- 
pute the dipole moment and the populations of an exactly 
solvable one-dimensional model. In Section III we map 
the interacting problem onto a non-interacting Kohn- 
Sham one and carry out numerical time propagations us- 
ing two different adiabatic approximations, ALDA and 
EXX. We also derive analytic expressions for the dipole 
moment and the populations in the KS system within 
the adiabatic EXX approximation and discuss and quan- 
tify the appearing dynamical detuning. We conclude our 
paper with a summary and outlook in Section IV. 



II. RABI DYNAMICS IN INTERACTING 
MANY-BODY SYSTEMS 

The Hamiltonian for an interacting non-relativistic A^- 
electron system coupled to a laser field of the form £{t) = 
£q sin(ajt) in dipole approximation is given by 



N 



£{t) 



(1) 



with T being the kinetic energy, Vext describing the exter- 
nal potential, and T4e being the electron-electron inter- 
action (atomic units m = e = Ti = 1 are used throughout 



the paper) . We denote the eigenstates and eigenvalues of 
Hamiltonian (1) with ip^ and e^, respectively. In order 
for Rabi's solution to be valid (for a detailed discussion 
see e.g. Ref. 7), the system under study needs to be an 
effective two-level system reducing the solution space to 
tpg (ground state) and V'e (dipole allowed excited state) 
with eigcncnergics eg and eg. The system then has a reso- 
nance at A — eg — eg. For a generic situation the two- level 
approximation should be valid if the laser frequency w is 
close to the resonance, and the amplitude of the driving 
field is not too strong to disturb the rest of the spectrum. 
These conditions can be formalized as follows 

5 « A, no«uj, (2) 

where (5 = — A is the detuning from the two-level res- 
onance and r2o = d^gSg is the resonant Rabi frequency, 
with deg = {tpg I I V'e) being the dipole transition ma- 

trix element. 

For an effective two-level system the time-dependent 
many-body state |^'(i)) can be written as a linear com- 
bination of both ground and excited states, i.e. 

|^(i)) =a,(i)|V'<,)+ae(i)l^e) (3) 

with |ag(t)p = ng{t) and |ae(t)P = ne{t) being the time- 
dependent populations of the ground and excited states, 
respectively. Normalization of the wave functions implies 
ng{t) + ne{t) = 1. Using Eq. (3) to compute the dipole 
moment d{t) ^ (VX^)! fjlV^t)) = deg2Re[a*g{t)ae{t)] 
we obtain 



d{t) = 2de,g^ng{t)ne{t)coii{ujt + ip{t)), (4) 

where ip{t) is a slowly varying (on the scale of 1/w) part 
of the phase difference of the coefficients ag{t) and ae(t). 
The Rabi dynamics is encoded in the envelope of the 
dipole moment through the time dependence of ng{t) and 
ne{t). After the two- level projection the time-dependent 
Schrodinger equation idt\ip{t)) = H\ilj{t)) reduces to the 
form 




The conditions (2) imply the existence of two well- 
separated timescales: one governed by the external fre- 
quency w and the other determined by the Rabi frequency 
fioj the frequency of the oscillations between the ground 
and the excited state. This timescale separation allows 
for the use of the rotating wave approximation (RWA) 
[7] leading to the following differential equation for (t) 
(for an outline of the derivation sec Appendix A) 

d^n,it) = ~{d^ + nl)n,{t) + ^nl (6) 

with initial conditions ^^(O) — and 9(71^(0) = 0. Eq. 
(6) describes a harmonic oscillator with a restoring force 
which increases with increasing detuning 6 (see potentials 
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FIG. 1. Dipole moment (red) and populations Tie (solid black 
line) and Ug (dashed black line) from the analytic solution of 
(5) using So = 0M25LU for detuning S = OMflo (0.0006 Ha) 
(a) and 5 = 2.2^0 (0.016 Ha) (b). 

for different detunigs S in Fig. 3). Thus, increasing 6 
results in a squeezing of the harmonic potential leading to 
a decrease of the amplitude of the oscillations according 
to n'l^^^ = r2Q/(ilQ + (5^) and a larger Rabi frequency (see 
Fig. I). 

In order to investigate the description of Rabi oscilla- 
tions within TDDFT wc analyze a one-dimensional (ID) 
two-electron model system which has the advantage that 
it can be solved exactly [8]. For ease of comparison we 
choose the same model as in [3] with the external poten- 
tial 

Ve^t = + l (7) 

and the electron-electron i nteraction bein g of soft- 
Coulomb type, i.e. T4e = 1/ a/ {xi — X2)^ + 1. The eigen- 
functions and eigenvalues are obtained by diagonalization 
of the Hamiltonian (1) using the octopus code [9, 10]. 
The calculations are performed in a box from —100 to 
100 bohr with a spacing of 0.2 bohr. The obtained eigen- 
values arc eg = —2.238 Ha and eg ~ —1.705 Ha, and 
the static dipole matrix element is deg = 1.104. In or- 
der to induce Rabi oscillations, a laser field of the form 
£{t) = So sin(wi) is turned on at t = 0, with Eq = 0.0125a; 
which ensures conditions (2) are satisfied. The frequency 
uj of the applied field has been chosen to be close to the 
resonance A = — eg. In Fig. 1 the time-dependent 
dipole moment d{t) and the populations ng(t) and Ugit) 
for 6 = 0.08f7o and S = 2.2^lo arc shown. The effect of 
the detuning manifests in an incomplete population of the 
excited state and a consequent decrease in the amplitude 
of the envelope of d(t) that is proportional to ^Jn^fQ. For 
small detuning the minima and the maxima of rtg coin- 
cide with minima of the envelope, but for larger detuning 
the dipole moment only goes to zero for the minima of 
rig. In Fig. la, 5 is very small but non-zero leading to the 
appearance of a neck at the odd minima of the envelope 
function which coincide with the minima of Ug. The neck 



increases with increasing 5 and evolves into a maximum 
for Fig. lb. Thus, the first minimum of Fig. lb corre- 
sponds to one complete cycle and can be identified with 
the second minimum in la. We note that looking only 
at the dipole moment is insufficient to discern between 
resonant and detuned Rabi oscillations, only studying 
the population ng of the excited state provides access to 
the complete picture. A comparison between the ana- 
lytic solution of Eqs. (4)-(6) and the results of the time- 
propagation with the octopus code shows perfect agree- 
ment, which confirms that the conditions (2) are fulfilled 
for the chosen values of £q and w. 

III. RABI OSCILLATIONS IN THE KOHN 
SHAM SYSTEM 

In TDDFT the interacting system is mapped onto a 
non-interacting KS system which reproduced the correct 
density p{r,t) [1, 2]. The time-dependent KS Hamilto- 
nian corresponding to (1) is given as 

Hs^H'J + V,^y:it)+r£{t), (8) 

where the static KS Hamiltonian reads H^! = T + Vext + 
14i2;c[po]- The KS wave functions 0j(r) are eigenfunc- 
tions of with eigenvalues and the time-dependent 
density is computed as p{r,t) = |(/)j(r,i)p. Here, we 
are studying a two-electron singlet system, but within 
the two-level approximation there is always one unique 
orbital that is getting deoccupied and another that is get- 
ting populated, independent of the number of particles. 
Thus, the time evolution only affects one unique orbital 
and follows from the KS equation idt(j>{r,t) = Hs(j){r,t) 
with initial condition (/)(r, t = 0) = 0g(r). The KS 
equation is non-linear due to the dependence of the 
Hartree-exchange-correlation potential Vhxc on the den- 
sity which, for two electrons in a singlet state, is given 
as p{r,t) = 2|(/)(r, The time-dependent dipole mo- 
ment d{t) is an explicit functional of the time-dependent 
density, i.e. d{t) = J cfirp{r, t)r. The exact KS system re- 
produces the exact many-body density p{r, t) and, hence, 
the exact dipole moment d{t). However, this need not 
be true for an approximate functional. Especially, us- 
ing adiabatic approximations has been shown to have a 
dramatic effect on the calculated density during Rabi os- 
cillations [3]. 

Here, we employ two different approximations to the 
one-dimensional xc potential, the recently derived adi- 
abatic local density approximation (ALDA) [8] and ex- 
act exchange (EXX) which for a two-electron singlet is 
adiabatic and equal to Hartree-Fock. The resonant fre- 
quencies are calculated from linear response in octopus 
which yields uj^^da ^ q ^ia and ljj'^^^ = 0.549 Ha. 
We then apply a laser field in analogy to the exact cal- 
culation with an amplitude of £q = 0.0125a; using the 
resonant frequency for each case. Propagating with the 
ALDA and EXX results in the dipole moments shown in 
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FIG. 2. Dipole moment (red) and populations (solid black 
line) and (dashed black line) for ALDA (a) and EXX (b). 
The results from the theoretical model Eq. (16) are given in 
(c). The calculations in (a) were performed for lu^'^^^ = 
0.476 Ha and the ones in (b) and (c) for oj'^^^ = 0.549 Ha. 

Fig. 2a and 2b respectively. Note that despite of the 
applied laser being in resonance with the system both 
EXX and ALDA approximations show the characteristic 
signature of detuned Rabi oscillations: the population of 
the excited state is incomplete and the zeros of the dipole 
moment coincide with the minima of but not with the 
minima of rt^, giving a similar picture as in Fig. lb. As 
the KS system only has to reproduce the exact density, 
we do not expect a complete population of the excited 
KS orbital. However, as we can sec from Fig. 2, the 
system remains mainly in its ground state and, like al- 
ready pointed out in [3] , its density resembles mostly the 
ground-state density. This does not imply that the KS 
system is not undergoing Rabi oscillations. Based on the 
behavior of the populations we postulate that the oscilla- 
tions seen in Fig. 2 are in fact detuned Rabi oscillations 
and not of classical origin as claimed in [3] . 

If an adiabatic approximation is used the potential at 
time t is a functional of the density at this time, i.e. 
yhxc{t) = yhxc[p{t)\- In the following, we show that 
Vhxc{t) indeed introduces a detuning that drives the sys- 
tem out of resonance. We again rely on the conditions 
(2). i.e. describe the KS system as an effective two- level 
system with 

0(i)(r,t)=a^(t)0,(r)-f a:(t)</)e(r). (9) 

Projecting the KS Hamiltonian (8) onto the two-level KS 
space (9) yields the 2x2 matrix 

/ el + ef{t) dl^8{t)+Fxc{t)\ 
\dl^S{t) + Fxc{t) e? + ef(t) ) ^'"^ 

with the dipole matrix element being 6?*^ — (0e|f|0g)- 



The additional terms, eg%{t) = {(f>g.e\V^xc i^)\'l>a-e) 

describe the fictitious time- 



it) = {<i^,Kt:me^ 

dependence that results in a dynamical detuning of the 
Rabi oscillations. As in the linear Rabi oscillations, the 
matrix (10) determines the coefhcients ag{t) and al{t) 
and the equations of motion for the dipole moment d'' (t) 
and the population Tig (<) = |ag(<)p. Compared to Eq. (5) 
we note that each entry in (10) contains an additional 
term depending on V/f^". In order to investigate the 
consequences of this term we again study the external 
potential (7) and use the EXX functional for which a rel- 
atively simple analytic expression for the additional ma- 
trix elements can be derived. The behavior using ALDA 
is very similar to the one for the EXX approximation (see 
Fig. 2a). The analysis, however, is more involved due to 
the functional not being linear in the density. 

For the two-electron singlet case investigated here, 
the Hartree-exchange-correlation potential V^^"^(a;, t) is 
equal to half the Hartree potential and, hence, given as 



hxc I d^xVeei\x-x'\){po{x')+Sp{x',t)). 



yEXX 



(11) 

Here, the part containing po determines Vh2:c[po] while 
Sp results in the additional V'^J'"- We then rewrite the 
contributions to the diagonal terms of Eq. (10) as 
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FIG. 3. Potentials corresponding to the diflerential equations 
(6) and (16). The dynamical detuning leads to a quartic po- 
tential which has a similar eflFect to the potential as a large 
detuning in the applied frequency (see text for details). 



differential equation for the level population nl{t) (see 
Appendix B for the derivation) 



nlit) + -ni (16) 



CW=^9<W. 6r(i) = Ae<(i) 



where, in EXX, the coefficient A reads 



(12) 



An.R — 



'g,e - I jdxdx' (\(l)e{x')\^ 



\M^')\')Vee{\x-x'\MgA^W 

(13) 

For the ofF-diagonal contributions we recall that d'^{t) ~ 
2dlgRe[ag{t)* al{t)], as in the exact case, and rewrite the 

contribution of V^^^{t) to the off-diagonal terms as a co- 
efficient g multiplied by the time-dependent dipole mo- 
ment 



d%t) 



"eg 



(14) 



Here, g is given as 
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dxdx (l)e{x')4>g{x)Vee{\x - X \)4>g{x)4>e{x) . (15) 



The coefficient g also enters when one calculates the res- 
onant frequencies in linear response. Within the two- 
level approximation the resonant frequency is given as 
Wq = \/ As (As + 2g) (sec Appendix B) which yields 
uj^^^ = 0.532 Ha. The deviation from uj^^^ calculated 
from time propagation of Hamiltonian (8) in octopus is 
of the order of 3%, coinciding with the deviation of our 
system from a true two-level system which we estimate 
from (l — (rigit) + nl(t))) . Using again the RWA we ob- 
tain, to leading order in X/loq and g/ujQ, the following 



with ils — dlgSo and 7 = A — 2^. Neglecting higher or- 
der terms leads to an error of about 10%. Within these 
error bars we can safely use Eq. (16) because it contains 
all the relevant physics. Unlike Eq. (6) which represents 
a harmonic oscillator, Eq. (16) corresponds to an anhar- 
monic quartic oscillator and its solution can be written 
in terms of Jacobi elliptic functions [11]. Equivalently, 
Eq. (16) can be integrated numerically. Comparing Eqs. 
(16) and (6) we conclude that the adiabatic approxima- 
tion introduces a time-dependent detuning proportional 
to -^jnl{t). The anharmonic oscillator no longer repre- 
sents a parabola but the dynamical detuning also results 
in an increase of the restoring force. This can be clearly 
seen in Fig. 3 where we plot the potential correspond- 
ing to the restoring force in Eqs. (6) and (16) for the 
detunings of Fig. 1 and the EXX calculation( Fig. 2c). 

Using the same ID model system Eq. (7) as before 
we apply a field of amplitude Eq = 0.0125w^"^'^ and 
calculate the dipole moment d^{t) and population nl{t) 
from Eq. (16). For this system we obtain for the bare 
KS eigenvalues = -0.750 Ha and = -0.257 Ha 
which yields As = 0.494 Ha. For the various matrix el- 
ements we obtain d^g = 0.897, g = 0.071, A = -0.125, 
and 7 ~ —0.268. The results are shown in Fig. 2c in 
comparison to the numerically exact time-propagation in 
octopus (Fig. 2b). The discrepancy between the numer- 
ical propagation and the analytical results is mainly due 
to the fact that we kept only the leading orders in A/wo 
and g/u!o in the derivation of Eq. (16). However, our 
simple model clearly captures the effect of the dynami- 
cal detuning present in all adiabatic functionals. In the 



5 



physical system the density changes dramatically during 
the transition and, in order to be able to describe res- 
onant Rabi oscillations, a functional has to keep track 
of these changes over time, i.e. it needs to have memory. 
For any adiabatic functional the potential will change due 
to the changing density and the system is driven out of 
resonance. We emphasize that this effect is not limited 
to TDDFT but is generic for all mean-field theories, e.g. 
HP or all hybrids, when the effective potential depends 
instantaneously on the state of the system. 



IV. CONCLUSION 

In conclusion, we demonstrate that the use of adia- 
batic approximations leads to a dynamical detuning in 
the description of Rabi oscillations. Only the inclusion 
of an appropriate memory dependence can correct the 
fictitious time-dependence of the resonant frequency. In 
particular, by clearly identifying the reason behind the 
failure of adiabatic approximations, we are able to prove 
the corrcspronding conjecture recently made in [3]. Our 
results constitute a very stringent test for the develop- 
ment of new xc functionals beyond the linear regime as 
all (adiabatic) functionals available till now fail to re- 
produce Rabi dynamics. Adiabatic functionals will fail 
similarly in the description of all processes involving a 
change in the population of states. Work along the lines 
of deriving a new memory-dependent functional is ongo- 
ing. The description of photo-induced processes in chem- 
istry, physics, and biology and the new field of attosec- 
ond electron dynamics and high-intense lasers all demand 
fundamental functional developments going beyond the 
adiabatic approximation. 



ACKNOWLEDGMENTS 

We acknowledge support by MICINN (FIS2010-21282- 
C02-01), ACI-promociona (ACI2009-1036), Grupos Con- 
soHdados UPV/EHU del Gobierno Vasco (IT-319-07), 
the European Community through e-I3 ETSF project 
(Contract No. 211956), and the European Research 
Council Advanced Grant DYNamo (ERC-2010-AdG - 
Proposal No. 267374). JIF acknowledges support from 
an PPI-fellowship (FIS2007-65702-C02-01). 



Appendix A: Equation of motion for ne(t) in the 
interacting many-body problem 

We provide here a derivation of the equation of motion, 
Eq. (6), for the population ne{t) of the excited state. 

As a first step we formulate a system of equa- 
tions for physical observables, the dipolc moment 
d{t) = 2degRe[a*ae\, the "transition current" J{t) ^ 
2deglm[a*ae], and the populations Ug^eif) = |ag,eP- Us- 
ing the Schrodingcr equation Eq. (5) for ag,e{t) and the 



above definitions of d{t), J{t) and Ug^eit) one can de- 
rive the following coupled differential equations for these 
quantities, 

dtdit) = AJ(t), (Al) 
dtJit) = -Adit) - 2ing - n,)\d,g\''£{t), (A2) 
dtueit) - -J{t)£it), (A3) 



where A = ep 



€g with €g^e 



degSit) = {^e\H{t)\^g 



{ipg,e\H{t)\ipg^e) and 
Equations (Al), (A2), and 
(A3) correspond to the two-level version of the conti- 
nuity equation, the equation of motion for the transition 
current, and the energy balance equation, respectively. 
These equation have to be supplemented with the nor- 
malization condition ng(t) + ne(t) = 1, and the initial 
conditions d(0) = J(0) = ne(0) = 0. 

By combining Eqs. (Al) and (A2) one obtains an equa- 
tion of motion for the dipole moment d{t) 

dfd{t) = -A^d{t) - 2(1 - 2ne)A\deg\^S{t). (A4) 

The next step in the derivation is to simplify this equa- 
tion using RWA. Namely, we separate "fast" and "slow" 
time scales by writing the dipole moment in the following 
form 



d{t) = deg[bi{t) cos{ujt) + 62 (t) sin(a;i)] 



(A5) 



where the coefficients bi{t) and 62(0 are assumed to be 
slowly varying {dtb/b ^ uj = A + d). Inserting the ansatz 
of Eq. (A5) into Eq. (A4) and making use of the condi- 
tion 6,Qa << A we arrive at the system of first order 
differential equations for slow variables 

dMt) - Sbi{t) = (A6) 
9t5i(t) + (562(i)-f^o(l-2ne) = (A7) 

From Eqs. (A3) and (Al) we obtain the equation of mo- 
tion for the populations 



dtUe = ioofei(t). 



(A8) 



where we again kept only leading in uj terms and, as usual 
in RWA, neglected irrelevant terms oscillating with 2lu. 
Combining this equation with Eq. (A6) results in 



dt 



, / N 2(5 
b2it)-—n,it) 

"0 



= 



which, together with the initial conditions ne(0) 
dtd{0) = 0, yields 

2S 

62 (t) = T^neit) 

"0 



and, due to Eq. (A6), 

biit) 



-dtUeit) 



(A9) 
and 

(AlO) 
(All) 



Finally, by inserting Eq. (AlO) and Eq. (All) into 
Eq. (A7) we arrive at the differential equation Eq. (6) 
for the population n,, (t) of the excited state in the many- 
body interacting system. 
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Appendix B: Equation of motion for nl{t) in the KS 

system 



Following the same scheme as in appendix A we do now 
derive Eq. (16) describing dynamics of the KS population 
nl{t). The equations of motion for the KS quantities, 
d^{t) = 2dlgRe[{airal], J%t) = 2d^/m[(a^)*<], and 
rig g(t) = jflg are similar to Eqs. (A1)-(A3) and can be 
obtained with the following replacements:A — > A(t) = 
A« + ef (i)-e^=(t), anddegfW ^ d^^S (t) + J^.,,{t) . One 
can again derive an equation of motion for the dipolc 
moment (in analogy to Eq. (A4)) which, due to the ad- 
ditional time dependence in A(t), aquircs an additional 
term and reads 



dtAjt) 

m 



dtd'it) 



(Bl) 



2(n^,-<)A(i)4g(d:/(t)+-F.c(0), 



where J-xc{t) is given by Eq. (14). For the two-electron 
singlet case and EXX approximation, the timc-depcndcnt 
contribution to the Hartrcc-cxchangc-corrclation po- 
tential is proportional to Sp{x,t). As 5p{x,t) = 
2nl{t){\(j)e{x)\^ - |0g(a;)p) + 2d"(i)0g(a;)(/)e(a;) we obtain 



Xnlit) 



(B2) 



with A = Ae — Ag and Ag e defined in Eq. (13). 

The resonant frequency uj^-^^ is now not given by the 
KS energy difference Ag but needs to be calculated from 
the linear response. In the linear regime Eq. (Bl) reduces 



to the following form 

d^,d^{t) = -A, ((A, + 2g)d^{t) - 2/^Mlg?m) , (B3) 

From the first term in the right hand side in Eq. (B3) we 
identify the linear response resonant frequency as 



.EXX 



= V(As+2g)A, 



(B4) 



We can now apply the same procedure to Eq. (Bl) as 
in the interacting case to derive an equation of motion 
for the occupation nl{t). Below, as well as in the main 
text, for the KS system we consider only the case of a 
resonant excitation, ur = uiq-^-^ , i. e. 6 = 0. Employing 
the ansatz of Eq. (A5) for the KS dipolc moment #(t), 
in analogy to Eqs. (A6) and (A7) we obtain 

dtbUt) + (A - 2.9) nl{t)b-lit) = (B5) 
dtbl{t)-{X-2g)nl{t)b-;{t) 

-as(l - 2<(t)) = (B6) 

with rts = dlgSo. From the EXX analog of Eq. (A3) 
(the energy balance equation) we find that to the leading 
order in RWA the equation of motion for the occupation 
has exactly the same form as Eq. (A8), i. e.. 



Therefore, with the initial conditions nl{0) 
dtnl{0) = we get for the two coefficients 



blit) 



A-2g 



Kit)? 



(B7) 
and 



(B8) 
(B9) 



Finally, by inserting Eq. (B8) and Eq. (B9) into Eq. (B6) 
we arrive at the differential equation Eq. (16) for the 
population ^^(t) of the KS excited state. 
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